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GALACTIC COORDINATES 


DECUS Program Library Write-up 


DECUS NO. 8-501 


ABSTRACT 

This program consists of two parts: the first converts star positions from astronomical 
coordinates (right ascension and declination) to galactic coordinates (latitude and longitude) 
the other converts in the reverse direction. 

USAGE 

To convert from astronomical coordinates (right ascension and declination) to galactic 
latitude and longitude, we will start with the following set of axes: 

+X1 AXIS THRU RA=6 HR 49 MIN, DEC=0 DEG 
+Y1 AXIS THRU RA=12 HR 49 Ml IN, DEC=0 DEG 
+Z1 AXIS THRU DEC=+90 DEG 


Then 


XI = COS DEC DOS (RA-6H49M) 

Y1 = COS DEC SIN (RA-6H49M) 

Z1 = SIN DEC 

We will now rotate these coordinates about the +X1 axis so that the Z2 axis is thru RA=12 
HR 49 MIN, DEC=27.4 DEG. The angle of rotation is Rl=62.6 DEG: 
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Now we have: 

+X2 AXIS THRU GALACTIC EQUATOR AT LOnGITUDE=21 3 DEG 
+Y2 AXIS THRU GALACTIC EQUATOR AT LOnGITUDE=3 < 03 DEG 
+Z2 AXIS THRU NORTH GALACTIC POLE 

Finally, we rotate thru an angle of 213 degrees, then calculate the galactic latitude and 
longitude. 

The program should be compiled, loaded and started in the usual manner. It will type: 

HRS & MINS OF RIGHT ASCENSION= 

and wait for the operator to key in these values with suitable terminators (spaces). Similarly 
the declination is entered. After a quarter second of calculation, the machine types out 
the answers and starts a new problem. 
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Ci FROM ASTRONOMICAL TO GALACTIC COORDINATES, BY A. 
a; FORMAT C / ,"HRS ft MINS OF RICHT ASCENSION*") 

3; FORMAT (E,E) 

a; FORMAT C/,"l)ECS ft MINS OF DECLINATION*") 

5; FORMAT C/,"GALACTIC LATITUDE=",I,"DECS ",I,"MINS’ 
"LONG ITUDE=",I,"DECS ", I,"MINS",/,/,/) 

RD=57.29578 
DR=.I745329E-1 
R=62.6*DR 
5R=SINFCR) 

CR=COSFCR) 

1 J TYPE a 

ACCEPT 3,D,EM 

H=C360.-15.*CD-6.)-.25*CEM-49.)>*0R 
TYPE 4 

ACCEPT 3,D,EM 
D= C D+EM/60•)*DR 
CD=COSFCD) 

21=SINFCD) 

Y1=CD*SINFCK) 

X2=CD*C0SFCB) 

Y2=71*CR+2l*SR 
22=-Yl*SR+2l*CR 
IF CX2) 7,8,9 

7I CL0N=33.-RD*ATNFCY2/X2) 

CO TO 15 

8; IF CY2) 11,12,13 
11 i CLON=123. 

CO TO 15 
12;IC=90.*22 
IH=MIN=M=3 
CO TO 14 
13;CLON=303. 

GO TO 15 

9; GL0N=213.-RD*ATNFCY2/X2) 

15;IF C CLON) 16,17,17 
16;CLON=GLON+360. 

17;IH=CLON 
D= IH 


MOSES, 20 OCT 71 


,/,"GALACTIC 


MIN=CCLON-D)*60. 

CL4T=8D*ATNFC Z2/3QTFCX2*X2+Y2*Y2)) 
I C = CLAT 
D= I C 

M=C CLAT-D)*60• 

14,‘TYPE 5, IC,M, IH,MIN 
CO TO 1 
END 
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c; FROM GALACTIC TO ASTRONOMICAL COORDINATES, BY A. MOSES, 20 OCT 71 
2J FORMAT (/, '’DECS & MINS OF CALACTIC LATITUDE 3 ") 

3J FORMAT CE, E) 

4; FORMAT (/, '’DECS & MINS OF CALACTIC LONGITUDE 3 ") 

5; FORMAT (/,“RIGHT ASCENSION 3 ",I,"HRS ",I,"MINS ", I,"TENTHS",/, * 
"DECLINATION 3 ",I,"DECS ", I,"MINS",/,/,/) 

RD=57.29578 
DR=.1745329E-1 
R=62.6*DR 
SR=SINFCR) 

CR=COSFCR) 

1; TYPE 2 

ACCEPT 3,D,EM 
D=C D+EM/60•)*DR 
CD=COSFC D) 

Z2=SINFCD) 

TYPE 4 

ACCEPT 3,D,EM 

D=CD+EM/60.>213.)*DR 

Y2=CD*SINFCD) 

Xt=CD*COSFCD) 

Y1=Y2*CR+Z2*SR 
Z1=-Y2*SR+32*CR 
IF CXI) 7,8,9 

7) H=18.81667-RD+ATNF(Y1/-XI)/!5. 

CO TO 15 

8J IF CY1) 11,12,13 
ll;H=12.81667 
CO TO 15 
12;IC=90.*Z1 

IH=MIN=M10=M=0 
GO TO 14 
13)H=.8166667 
CO TO 15 

9J H=6.816667+RD+ATNFC Yl/Xl)/15. 

15)IF CH-24.) 16,17,17 
17)H=H-24• 

16;IH=H 
D=IH 

EM= C H-D)*60• 

MIN=EM 

D=MIN 

Ml 0 = CEM-D)*10. 

D=RD*ATNFC Z1/SQTFCX1*X1+Y1*Y1)) 

I C = D 
EM= I C 

M=C D-EM)*60 * 

14;TYPE 5,IH,MIN,MI0,IC,M 
GO TO 1 
END 
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